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ABSTRACT: We study a modified two-dimensional dilaton gravity theory which is exactly 
solvable in the semiclassical approximation including back-reaction. The vacuum solutions 
of this modified theory are asymptotically flat static space-times. Infalling matter forms a 
black hole if its energy is above a certain threshold. The black hole singularity is initially 
hidden behind a timelike apparent horizon. As the black hole evaporates by emitting 
Hawking radiation, the singularity meets the shrinking horizon in finite retarded time to 
become naked. A natural boundary condition exists at the naked singularity such that 
for general infalling matter-configurations the evaporating black hole geometries can be 
matched continuously to a unique static end-state geometry. This end-state geometry is 
asymptotically flat at its right spatial infinity, while its left spatial infinity is a semi-infinite 
throat extending into the strong coupling region. 
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Hawking's discovery that black holes radiate thermally gave rise to a long- 

standing question concerning the consequences of combining quantum theory and general 
relativity. [@,@,@,0] Does evolution from an initial pure state take place unitarily to a 
final pure state or non-unitarily to a final mixed state? Intimately linked to this question 
is the final geometry resulting from black hole evaporation. 

Here we present a specific two-dimensional (2D) dilaton gravity model in which a 
black hole evaporates leaving a static semi-infinite throat as the end-state or "remnant" 
geometry. Our model is a modification of the CGHS model. [||] We solve the semiclassical 
equations and get closed-form expressions for the metric and dilaton field. 

The classical 2D CGHS action H is 



Sa = ^ d 2 xV^ (R^ + 4(V</>) 2 + 4A 2 ) - \ £(V/ 

71 J L i=l 



where </> is the dilaton field, is the 2D Ricci scalar, A is a positive constant, V is 
the covariant derivative, and the fi are iV matter (massless scalar) fields. The action (]!]) 
describes a 2D effective theory in the throat region of a 4D almost extreme magnetically 
charged black hole.[[]9[],[]10[]] It may also be regarded as a 2D arena in which some of the 
main questions about black hole evaporation can be studied. Among the classical solu- 
tions stemming from the action (P are vacuum solutions, static black hole solutions, and 
dynamical solutions describing the formation of a black hole by collapsing matter fields. 
For a review see Ref. 

To study one loop quantum corrections and back-reaction one can use the trace- 
anomaly for massless scalar fields in two dimensions, (T^) = -J^R( 2 \ and find the effective 



action Spl for which (T MI/ ) = — -^L^j^SpL- This is the Polyakov-Liouville action [||T2 



Spl = J d 2 xy^(x) J d 2 x' v ^{x 1 )R i2 \x)G(x,x')R {2 \x'), (2) 

where G(x,x') is a Green function for V 2 . Here we take the large iV limit, in which 
h goes to zero while N% is held fixed. In that limit the quantum corrections for the 
gravitational and dilaton fields are negligible, and one need take into account only the 
quantum corrections for the matter (scalar) fields. The one-loop effective action is then 
Si\\ = S c £ + NSpL. There are no known analytic solutions to this one-loop effective theory, 
though there are some numerical ones.[|13[]] In order to find analytic solutions including 
semiclassical corrections, one can modify the action as in [ ||14|| , ||15|l )|1l6n. Our approach 
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is similar, in that we modify the original CGHS action ([!]) and find analytic solutions 
to the modified equations including back-reaction. However, our analytic solutions yield 
closed-form expressions for the metric and dilaton field. This allows us to fully analyze the 
solutions. 

We add to the classical action ([1]) a local covariant term of one-loop order, 



Nfi 

corr 24tt 



j <fV=0 ((V0 2 - #i< a) ) . (3) 
Now the total modified action including the one-loop Polyakov-Liouville term is 

'-'mod = c ^ ^ '-'corr 

+ NS PL . (4) 

Using null coordinates z ± and conformal gauge g ++ = = 0, = — \e 2p (ds 2 = 

—e 2p dz + dz~), the action can be written in the form 

1 r r ~\ N 

S mod = - dz+dz ~ 2d-(0-p)d + (e- 2 *-^(cf ) -p))+\ 2 e 2 ^ + -J2d + fid-f l 

(5) 

where k = [From the point of view of string theory, the action (^) with free fields X = 
e~ 2 ^ and Y = <ft — p (which are flat target space coordinates), describes a conformal field 
theory with tachyon and dilaton backgrounds T = — 4A 2 e _2y and $ = — 2X + 2kY.\^ t ?]]] 



The action (|^) is also invariant under the transformation^ |fT4"f 5(p = 5p = ee 2< ^, with the 
conservation equation <9 M <9^(</> — p) = 0. We therefore can complete the gauge fixing by 
choosing the "Kruskal coordinates," x ± (z ± ), in which (fi(x + ,x~) = p{x + ,x~). In this 
Kruskal gauge the equations of motion derived from the modified action (|5|) are exactly 
the same as the classical ones 

9 x+ d x - ( e -2p(- + ,--)) = Q x+Q ^ ^-2^(x+,x-) j = _ A 2 (g) 

d x+ d x ~h{x + ,x~) =0, (7) 

while the constraints get modified by non-local terms t±(x ) arising from the Polyakov- 
Liouville action. In conformal gauge, one can use the trace anomaly of iV massless scalar 



t Unlike in the RST model, [14|, in this model the transformation is exactly the same as in the 
classical case. 
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fields fi to obtain (T^__) = —nd + d-p and integrate [[IB|,[I^|,[SDH the equation V At (T / { !/ ) = 



to get the quantum corrections to the energy-momentum tensor of the fi matter fields 

(T£ ± ) = K {d 2 ± p - (d ±P f - t ± {z±)) , (8) 

where t±(z ) are integration functions determined by the specific quantum state corre- 
sponding to the expectation value (^\Tf u \^) = (Tf). These functions can be determined 
by boundary conditions. Alternatively, Eq. (||) can be obtained by varying NSpl. Then 
the functions t±(z^) arise from the homogeneous part of the Green function in Eq. (||). 
Our modified constraints (in Kruskal gauge) are 



6S. 



mod 



±± 



o 



-dl ± I e 



-%4>{X+ ,X-)\ _ (rpf 



(T±±) c e + Kt±(x ) = 0, 



(9) 



where (T±±) c i = ^J2iLi(^x ± fi) 2 is the classical (zero order in K) contribution to the 
energy-momentum tensor of the fi matter fields. (T^) in (|8]) is the one- loop quantum 
correction of order h, so the full energy- momentum tensor of the /-fields is (7/f l/ ) c £ + 
(Tp + 0(h 2 ). 

For a given classical matter distribution and a given t±(x ) one finds the solution for 
the equations of motion (|6|) with the constraints (0): 



-2p 



-X 2 x + x — I dxt I dx± (T^ + ) c e — Kt + (x±) 



(10) 



dxl 



dx x (Ti_) c £ - Kt-{x x ) 



+ a+x + + d-X + b 



where a± and b are constants. First, let us consider the linear dilaton flat space-time 
solution, e -2 ^ = e~ 2p = — X 2 x + x~ . It corresponds to the choice (T^) c £ = and t±(x ) = 
a± = b = 0. To determine the corresponding quantum state one must calculate 
(t£ ± ) in (H) using the given t±(x ± ). In flat coordinates a^, which are related to the 
Kruskal coordinates x^ 1 by the conformal coordinate transformation ±Ax^ 



e ±Xa± , the 



We see that unlike in the RST model, in our 



expectation values (§) are (T£ ± (cr ± )) 

model (T±±(a :iz )) 7^ for the linear dilaton solution. Because (2± ± ) = and = 0, 

the quantum state corresponding to the linear dilaton solution may describe a system 
in thermal equilibrium at temperature T = 

In our model we also have static black hole solutions. [^T|] These correspond in Eq. (|T(i| ) 
to the choice (T^) c i = t±(x ± ) = a± = and b = M/X. For these solutions at future and 



past null infinity (3 + and respectively) one has (T± ± ) = ^j-; the solutions evidently 
describe a black hole in thermal equilibrium at temperature^ = ^- . This is as we would 
expect: A static black hole solution in a self-consistent semiclassical theory of Hawking 
radiation including back-reaction is possible only if the black hole is in thermal equilibrium 
with incoming radiation. 

In order to find the solution corresponding asymptotically to the Minkowski vacuum 
we can use (H) to find the solution for which (t£ ± (ct ± )) = 0. The functions t±(x ± ) 
are determined by imposing appropriate boundary conditions on 3^. We assume that 
on these boundaries the metric is flat, such that p(cr ± ) and its derivatives vanish in the 
asymptotically flat coordinates a . Then the first two terms on the right- hand-side of (||) 
vanish on the boundary and we get 

(^^boundary = 

We see from (^Tj) that the Minkowski vacuum corresponds to t±(cr ± ) = 0. To find the 
corresponding t±(x ) in "Kruskal coordinates," one can use the tensor transformation of 
(t£ ± ) in Eq. (under a conformal coordinate transformation) and get 

^) = (^)'(^ (0±) -5^^) = ^ (12) 

where Dy[z] is the Schwarz operator Dy[z] = dyZ/(d y z) — ^(d 2 z/d y z) 2 and we use 
t±(a ± ) = 0. Using (|T0|), fll2|) and (T^ u ) c £ = 0, we find that the general asymptotically 
Minkowski vacuum solution is 

e -20 = e -2 P = _ A 2 X + X - _ ^i g(_A 2 x+x-) + C, (13) 

where C is a constant. In asymptotically flat coordinates cr^ = t ± a, we have 

ds 2 = (1 - e~ 2Xa {K\a/2 - C))~\-dt 2 + da 2 ) 
<K<t) = -Xa + log (l - e- 2 ^(^a - C)) . ^ 

+ Since in 2D the Hawking temperature is mass independent, one may regard the linear dilaton 
solution as the zero mass limit of the static black hole solutions. This may explain the non-zero 
temperature of the linear dilaton solution in our model. 
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This solution is static, depending on the spatial coordinate a alone. On the boundaries 
3 s , the solution approaches the linear dilaton flat space-time solution, justifying our as- 
sumption. The reason this solution with no radiation at 3^ and the earlier ones with 
radiation there all asymptotically approach the linear dilaton flat space-time solution is 
that the coupling, e 2< ^, of the matter to the geometry vanishes exponentially fast at Of. 

Before we turn to the question of the ground-state solution, let us consider the ADM 
masses of the various solutions we have found. Suppose that we can choose as our ground- 
state one of the radiationless solutions fll4| ) with C = Co, where Co is a constant yet to be 
determined. Then the ADM mass [[22|, ||23|| ] of any other static solution ( |14"D is A(C — Co). 
On the other hand, the ADM mass of the linear dilaton solution as well as the static black 
hole solutions (relative to this ground-state) is infinite. This is already clear from the fact 
that these solutions have non-vanishing radiation on and can be checked explicitly by 



using the ADM mass definition. ||22j| These considerations make it plausible that one of 
the static solutions ([14]) should be the ground-state. We will see later that there exists a 
natural lower limit on Co which gives the preferred ground-state of lowest ADM mass. 

We next turn to the dynamical scenario in which the space-time is initially described 
by one of the static solutions in fll4]) (not necessarily the ground state solution Co), and 
in which a black hole is formed by collapsing matter fields. First we consider the simple 
shock wave solution, but our results can be easily extended to general infalling matter 
configurations. The shock wave of infalling matter is described by (t{_,) c £ = -^ T 5(x + — Xq) 
and (T[_) ce = 0.@] Unlike in the RST model, here we have a general initial static geometry, 
and the shock wave forms a black hole only if M, the energy of the shock wave, is above a 
certain threshold energy. We assume that M is above that threshold. Integrating (Tjf _i_ ) c £ 
in (|10|) and using (jl^) and a± = 0, we find the evaporating black hole solution 

k M 
e -2<t> = e -2 P = _ x z x + x - _ -log(-A 2 x + £-) T (x + - x+)0(x+ - x+) + C, (15) 

where O(x) is the standard step function. 

Before the shock wave, i.e., in the region x + < Xq, we have a static solution ([14]) which 
is not globally flat. If f [1 — log(fi:/4)] + C < 0, then the scalar curvature diverges on a 
timelike curve a = cr s , for which e - 2 ^ (<Tfl) = 0. Of course this is a region of strong coupling, 
and one would expect to have higher order quantum corrections there. On the other hand, if 
f [1— log(K/4)]+C > 0, the scalar curvature is bounded. Then the region {x + > 0, x~ < 0} 
is geodesically incomplete and one can analytically extend it to x~ > and x + < 0. Also 
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in this case there is a region of strong coupling near a = a m in = — 2^1°g(f )• I n the 
semiclassical approximation, one avoids the strong coupling region by imposing boundary 
conditions on a suitable time-like hypersurface. [[p4j] ,14, [^5| ] For the static solutions ( [14]) , 



(T^_ | _(a ± )) and (T/_(cr ± )) are constant on any time- like hypersurface a = const. Moreover 
(I+ + (c7 ± ))| ct=cto = (TLi^h^o fOT an y constant a . (16) 

This means that we can limit our model to a region in which the semiclassical approxima- 
tion is valid by imposing reflecting boundary conditions (|16|) on any time-like hypersurface 
a = oo that lies outside the region of strong coupling (these boundary conditions are 
also conformal [[25)). The geometry before the shock wave is therefore a static geometry, 
defined in the region a > a s in the case f [1 — log(K/4)] + C < (or a > a m i n in the 
case f [1 — log(ft/4)] + C > 0), with reflecting boundary conditions on a = a s + 5 (or on 
c = o'min + 5) where 5 is an arbitrary small positive constant. 

The solution to the future of the shock wave (x + > Xq ) is (see fll5D ) 

k M 

e -24> = e -2 P = _ A 2 X +( X - + A) - -log(-A 2 £ + £-) + — + C, (17) 

where A = ^+ . This solution is asymptotically flat and describes a black hole with a 
singularity at e~ 2c ^ = 0. The black hole singularity curve is 

k A/f 
-A 2 x+(xJ + A) - I log(-A 2 x^xJ) + — + C = 0. (18) 



Initially the singularity is behind an apparent horizon <9 + e 2< ^ = 0,[[26|,[10|] which is the 
curve 

-A a 4(afc+A) = £. (19) 

When the apparent horizon is formed, the black hole starts radiating. One can see this by 
calculating at future null infinity (x + — > oo). From ( |i~7D we see that the asymptot- 

ically flat coordinates on S + are a , related to cc ± by the conformal coordinate transfor- 
mation, Xa + = log(Xx + ) and — Xa~ = log(— X(x~ + A)). Using (JTJ) and ([12]) we get 

<^)>U = ^ (l - ■ (20) 



This is the "standard" Hawking radiation in 2D, where the Hawking temperature Tjj 
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it- is a constant. H One can further verify that when the black hole evaporates over a 



long period of time, i.e., if M >> kX, the spectrum of the Hawking radiation is indeed 
Planckian.[§,[0] 

As the black hole evaporates by emitting Hawking radiation, the apparent horizon 
shrinks and eventually meets the singularity in a finite proper time. They intersect at (see 
Fig. 1) 

1 / ( 4(M+>C) s K\ _ _ . / K ( 4(M + AC) n \ -1 

^t = Ti7 ^ +1) -t and x int = -A ( 1 - -e~ ( ^ +1) J . (21) 



A 2 A V 4 



At this point the singularity becomes naked. We show that it is possible to impose a 
boundary condition in which a weak shock wave emanates from the intersection point, 
resulting in a solution that is stable (having non-negative ADM mass), conserves energy, 
and is continuous with the metric defined to the past of the null hypersurface x~ = x~ nt . 

Before considering the solution to the future of the null hypersurface x~ = x~ nt (the 
end-state solution), we calculate the total amount E ra d of energy radiated during the 
evaporation. Integrating (|20|) over S + (up to x~ nt ) gives 

(TL(S-))da- = M + XC-^ (log(«/4) - 1) - (22) 

where a~ nt = a~(x~ nt ). The result (p2] ) is exact. The ADM mass [^2J of the dynamical 
solution ( ^5| ) (relative to the ground state C = C ) is Madm = M + X(C — C ). We 
see that the black hole radiates almost all of its initial energy. The unradiated mass 5M 
remaining as x~ -» x~ nt (which is the Bondi mass) is 

kX kAA 
5M = M ADM - E rad = — (log(«/4) - 1) - AC + ■—. (23) 

We now consider the solution to the future of the point of intersection (xf nt , x~ nt ). A 
natural candidate for such an end-state in our model is one of the static solutions ( [14]) , 
so we try to find boundary conditions such that the solution ([17]) is continuously matched 
to one of the static solutions (|14D • Remember that the asymptotically flat coordinates are 
a ± , so one should replace a in ([14]) with a = \(d + — a~). In the x ± coordinates the 
corresponding static solution is (see (|T^)) 

e -2<f> = e -2 P = _x 2 x+(x- + A) - ^log(-A 2 x+(x- + A)) + C. (24) 

We would like to see if there exists a constant C = C* , such that on the null hypersurface 
the solutions ([[7]) and ( |24"D can be matched continuously. This is indeed the 



int 
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case and from (|21]), (|l7|) and (f?§ we get C* = -~ (1 - log(«/4)). The end-state solution, 
or "remnant", is therefore 

e -24> = e -2 P = _ x z x +( x - + A) - ^log(-A 2 x+(x- + A)) - ^(1 - log(«/4)), (25) 

where x~ > x~ nt . From the constraint equations @ we find that 

1 N XA 
(TL(S-)U = p E( a -^) 2 = - Knt)- (26) 

This describes a shock wave originating at the intersection point and carrying a small 
amount of negative energy, K\A/(4x~ nt ), to null infinity. One may call it a "thunder- 
pop" . |I|J] The solution (^5|) is one of the static solutions that is asymptotically flat (with 
no radiation) on This means that there is no Hawking radiation after the thunderpop 

(ID- 

The mass remaining after the shockwave ( |26| ) has been emitted is 5M — k\A / \4x^ nt ) . 
One readily verifies that this is equal to the mass of the "remnant" (relative to Co) 
M rem = X(C* — Co). The fact that energy is exactly conserved, including terms of order 
h, supports the self-consistency of our semi-classical theory. Notice that C* and therefore 
the "remnant" mass is independent of the mass M of the infalling matter and of the con- 
stant C describing the initial static geometry. Even more surprising is the fact that the 
end-state solution with C = C* is the critical solution separating singular and non-singular 
static solutions described by Eq. (^3|). For C > C* the curvature of the solution ( pl|) is 
bounded, while for C < C* the curvature diverges on a time-like curve a = a s , for which 
e -24>(a B ) _ o j n so i u tion space, the solution ([25]) that has C = C* is the boundary between 
these two different classes of solutions. 

Consider the late-time space-like hypersurface £ shown in Fig. 1. Its right boundary 
(a — > oo) is z°, while its left boundary is the curve a = a cr , for which e~ 2< ^ = 0. For 
the critical solution we have d x +(e~ 2 ^ a ~ cr * > ) = e -2 ^ "^ = and the curve a = a cr is 
the analytical continuation of the apparent horizon to the region x~ > x~ nt . We define 
e = a — cr cr , and calculate the metric near e = 0. From (|25| ) we get 



2 -dt 2 + de 2 

where t = |(? + + &~)- An important feature of ( $7\) is that there is no linear term (in 
e) in its denominator. The first non- vanishing term is of order e 2 , which means that the 
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geometric structure near e = is that of an infinite throat. Consider for example the 
distance along t = constant curves. The distance to a = a cr diverges logarithmically, 
exactly as it does in higher-dimensional extremal black holes. The end-state space-time 
is geodesically complete. One may consider this solution as an "extremal 2D black hole" . 
On <7 = a cr the Ricci scalar is constant, = 4A 2 and the geometry is regular. 

The most natural choice of Cq for ground-state solution is the one with Cq = C* . 
This solution describes a static radiationless geometry which is regular everywhere. Any 
solution with smaller ADM mass (C < C*) has a naked singularity. In the class 
of solutions with no naked singularities, C = C* is the one with lowest energy. This is 
very similar to the linear dilaton vacuum solution (LDV) in classical dilaton gravity or to 
Minkowski space in Einstein gravity. § Also if we choose Co = C* , then the mass remaining 
after the thunderpop is exactly zero. Thus the end-state solution ([23]) is the static 
ground-state. Its geometrical structure is independent of the initial conditions and is a 
semi-infinite throat extending into the strong coupling region. 

In our 2D semiclassical model, one does not recover all the information of the initial 
state from the end-state solution. For infalling matter described by a general (r/ + ) c ^ of 
compact support, the solution ( JTH ) will depend only on the first two moments of (T^ + ) c £, 
M = Xf x + (Tl + ) cl dx+ and P + = J cl dx + .0 The end-state solution will still be 
(pop, but with A = X~ 2 P + . The information encoded in this "remnant" (or more precisely, 
in its past null boundary x~ = x~ nt ) is only about P + and M. Thus in our semiclassical 
model this end-state solution does not qualify as the "cornucopion" of Ref. [ PB|] ]. However, 
the semi-infinite throat extends to a region of very strong coupling. There may be sufficient 
freedom in this strong coupling region to encode more information through strong quantum 
gravitational effects. 

In this work we constructed an action in 2D dilaton gravity and showed that, with 
a natural boundary condition, all evaporating black holes in our model end in a unique 
ground-state geometry having a semi-infinite throat. 



Acknowledgements 

We thank Bruce Allen, John Friedman, and Jorma Louko, for helpful discussions and 
the National Science Foundation for support under grant PHY-9105935. 

§ Classical solutions with ADM mass smaller than the LDV have a naked singularity, as do 
Schwarzschild solutions with mass smaller than zero. 
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ABSTRACT: We study a modified two-dimensional dilaton gravity theory which is exactly 
solvable in the semiclassical approximation including back-reaction. The vacuum solutions 
of this modified theory are asymptotically flat static space-times. Infalling matter forms a 
black hole if its energy is above a certain threshold. The black hole singularity is initially 
hidden behind a timelike apparent horizon. As the black hole evaporates by emitting 
Hawking radiation, the singularity meets the shrinking horizon in finite retarded time to 
become naked. A natural boundary condition exists at the naked singularity such that 
for general infalling matter-configurations the evaporating black hole geometries can be 
matched continuously to a unique static end-state geometry. This end-state geometry is 
asymptotically flat at its right spatial infinity, while its left spatial infinity is a semi-infinite 
throat extending into the strong coupling region. 
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Hawking's discovery that black holes radiate thermally [1,2,3] gave rise to a long- 
standing question concerning the consequences of combining quantum theory and general 
relativity. [4,5,6,7] Does evolution from an initial pure state take place unitarily to a final 
pure state or non-unitarily to a final mixed state? Intimately linked to this question is the 
final geometry resulting from black hole evaporation. 

Here we present a specific two-dimensional (2D) dilaton gravity model in which a 
black hole evaporates leaving a static semi-infinite throat as the end-state or "remnant" 
geometry. Our model is a modification of the CGHS model. [8] We solve the semiclassical 
equations and get closed-form expressions for the metric and dilaton field. 

The classical 2D CGHS action [8] is 

1 r T 1 N 

Sd = ± d'x^-g (R^ + 4(V</>) 2 + 4A 2 ) - ± £(V/, 

71 J \_ i= i 

where cj> is the dilaton field, is the 2D Ricci scalar, A is a positive constant, V is 

the covariant derivative, and the fi are N matter (massless scalar) fields. The action (1) 
describes a 2D effective theory in the throat region of a 4D almost extreme magnetically 
charged black hole. [9, 10] It may also be regarded as a 2D arena in which some of the 
main questions about black hole evaporation can be studied. Among the classical solu- 
tions stemming from the action (1) are vacuum solutions, static black hole solutions, and 
dynamical solutions describing the formation of a black hole by collapsing matter fields. 
For a review see Ref. [11]. 

To study one loop quantum corrections and back-reaction one can use the trace- 
anomaly for massless scalar fields in two dimensions, (Tff) = ^i?*- 2 -*, and find the effective 

-g <V 



action Spi for which {T^) = — 7= c Spi- This is the Polyakov-Liouville action [12] 



S PL = J d 2 xy / -g(x) J d 2 x'y / -g(x')R (2 \x)G(x } x')R (2 \x') } (2) 

where G(x, x') is a Green function for V 2 . Here we take the large N limit, in which h goes to 
zero while NTi is held fixed. In that limit the quantum corrections for the gravitational and 
dilaton fields are negligible, and one need take into account only the quantum corrections 
for the matter (scalar) fields. The one-loop effective action is then = S c i + NSpl- 
There are no known analytic solutions to this one-loop effective theory, though there 
are some numerical ones. [13] In order to find analytic solutions including semiclassical 
corrections, one can modify the action as in [14,15,16]. Our approach is similar, in that we 
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modify the original CGHS action (1) and find analytic solutions to the modified equations 
including back-reaction. However, our analytic solutions yield closed-form expressions for 
the metric and dilaton field. This allows us to fully analyze the solutions. 

We add to the classical action (1) a local covariant term of one-loop order, 

Scorr = ^ / J 2 xV^ ((W) 2 - <f>R^ 
Now the total modified action including the one-loop Polyakov-Liouville term is 

'-'mod = Sci + 5"corr + NSpl- (4) 

Using null coordinates z^ 1 and conformal gauge g++ = g = 0, g+- = — \e 2p (ds 2 = 

— e 2p dz + dz~), the action (4) can be written in the form 

1 r r 1 N 

S mod = ~ dz+dz- 2d-^-p)d + {e- 2 ^-^-p))+\ 2 e 2 ^ + -Y J d + f l d-f l 
71 J \_ i= i 

(5) 

where n = -y^-. [From the point of view of string theory, the action (5) with free fields X = 
e~ 2 ^ and Y = cf> — p (which are flat target space coordinates), describes a conformal field 
theory with tachyon and dilaton backgrounds T = — 4A 2 e _2Y and $ = —2X + 2«y.[17]] 
The action (5) is also invariant under the transformationt[14] S(f> = Sp = ee 2 ^, with the 
conservation equation d^d^((j) — p) = 0. We therefore can complete the gauge fixing by 
choosing the "Kruskal coordinates," in which (f>(x + , x~) = p(x + , x~). In this 

Kruskal gauge the equations of motion derived from the modified action (5) are exactly 
the same as the classical ones 

d x+ d x - ( e -M* + ,*- A = dx+dx _ ( e -2<K* + ,*~)) = _A 2 (6) 



d x +d x -f t (x + ,x ) = 0, (7) 

while the constraints get modified by non-local terms t±(x^) arising from the Polyakov- 
Liouville action. In conformal gauge, one can use the trace anomaly of N massless scalar 
fields fi to obtain (T.£_) = —nd+d-p and integrate [18,19,20] the equation V' i (T / { J/ } = 
to get the quantum corrections to the energy-momentum tensor of the fi matter fields 

(Tl ± )= K (d 2 ± p-(d ± p) 2 -t ± (z±)), (8) 



t Unlike in the RST model, [14], in this model the transformation is exactly the same as in the 
classical case. 
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where t±(z^) are integration functions determined by the specific quantum state \^>) corre- 
sponding to the expectation value (\E r |T / { I/ |\E r ) = (T^ v ). These functions can be determined 
by boundary conditions. Alternatively, Eq. (8) can be obtained by varying NSpl- Then 
the functions t±(z^) arise from the homogeneous part of the Green function in Eq. (2). 
Our modified constraints (in Kruskal gauge) are 

^Sfi = -d 2 x± [e~^ + ^] - (Tl ± U + «<±(x ± ) = 0, (9) 

where (t£±) c i = ^ ^2iLi (d x ± fi) 2 is the classical (zero order in h) contribution to the 
energy-momentum tensor of the fi matter fields. (T^ v ) in (8) is the one-loop quantum 
correction of order h, so the full energy-momentum tensor of the /-fields is (T£ v ) c i + 
(Tp + 0(h 2 ). 

For a given classical matter distribution and a given t±(x^) one finds the solution for 
the equations of motion (6) with the constraints (9): 



e 



— c — /\ x x ( dx 2 i dx 



x 



dx r) / diX 



{TL_) c £ - Kt 



(10) 

+ a+x + + a-X~ + b 



where a± and b are constants. First, let us consider the linear dilaton flat space-time 
solution, e -2 ^ = e~ 2p = — X 2 x + x~ . It corresponds to the choice (Tj[ v ) c £ = and t±(x^) = 
a± = b = 0. To determine the corresponding quantum state one must calculate 
(t£±) in (8) using the given t±(x ± ). In flat coordinates cr 1 * 1 , which are related to the 
Kruskal coordinates x^ 1 by the conformal coordinate transformation ±Xx^ = e ±A(T± , the 
expectation values (8) are { ^ ( c" ^ ) ) = We see that unlike in the RST model, in our 

model (T£±(a :il )) ^ for the linear dilaton solution. Because (T£±) = and (T£_) = 0, 
the quantum state corresponding to the linear dilaton solution may describe a system 
in thermal equilibrium at temperature T = 

In our model we also have static black hole solutions. [21] These correspond in Eq. (10) 
to the choice (T£ v ) c i = t±(x^) = a± = and b = M/X. For these solutions at future and 
past null infinity and $5", respectively) one has (T£±) = the solutions evidently 
describe a black hole in thermal equilibrium at temperature+ T^h = This is as we would 



+ Since in 2D the Hawking temperature is mass independent, one may regard the linear dilaton 
solution as the zero mass limit of the static black hole solutions. This may explain the non-zero 
temperature of the linear dilaton solution in our model. 
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expect: A static black hole solution in a self-consistent semiclassical theory of Hawking 
radiation including back-reaction is possible only if the black hole is in thermal equilibrium 
with incoming radiation. 

In order to find the solution corresponding asymptotically to the Minkowski vacuum 
we can use (8) to find the solution for which ^±±((7^)) = 0. The functions t±(x^) 
are determined by imposing appropriate boundary conditions on Q^. We assume that 
on these boundaries the metric is flat, such that and its derivatives vanish in the 

asymptotically flat coordinates a^ 1 . Then the first two terms on the right-hand-side of (8) 
vanish on the boundary and we get 

(TL(^))\ hoimim , = (ID 

We see from (11) that the Minkowski vacuum corresponds to t±((7^) = 0. To find the 
corresponding t±(x^) in "Kruskal coordinates," one can use the tensor transformation of 
(T±±) in Eq. (8) (under a conformal coordinate transformation) and get 



'^' ± » = (|S) 2 (^ ± »-5 fl -^) = (SW- (12) 



where -D^[V] is the Schwarz operator Dy [z] = dyz/(d y z) — ^(d 2 z / 'd y z^ 2 and we use 
t±((7^) = 0. Using (10), (12) and (Tj[ v ) c £ = 0, we find that the general asymptotically 
Minkowski vacuum solution is 

e -24> = e -ip = -\? x + x - _ ^i og (_A 2 x+x-) + C, (13) 
where C is a constant. In asymptotically flat coordinates a^ 1 = t ± <r, we have 



ds 2 = (1 - e- 2X<7 (K\a/2 - C)) \-dt 2 + da 2 



4>(a) = -\a + log (l - e- 2 ^(^a - C)) . ^ 

This solution is static, depending on the spatial coordinate a alone. On the boundaries 
Q^, the solution approaches the linear dilaton flat space-time solution, justifying our as- 
sumption. The reason this solution with no radiation at and the earlier ones with 
radiation there all asymptotically approach the linear dilaton flat space-time solution is 
that the coupling, e 2 ^, of the matter to the geometry vanishes exponentially fast at Q^. 

4 



Before we turn to the question of the ground-state solution, let us consider the ADM 
masses of the various solutions we have found. Suppose that we can choose as our ground- 
state one of the radiationless solutions (14) with C = Co, where Co is a constant yet to be 
determined. Then the ADM mass [22,23] of any other static solution (14) is A(C — Co). 
On the other hand, the ADM mass of the linear dilaton solution as well as the static black 
hole solutions (relative to this ground-state) is infinite. This is already clear from the fact 
that these solutions have non- vanishing radiation on and can be checked explicitly by 
using the ADM mass definition. [22] These considerations make it plausible that one of 
the static solutions (14) should be the ground-state. We will see later that there exists a 
natural lower limit on Co which gives the preferred ground-state of lowest ADM mass. 

We next turn to the dynamical scenario in which the space-time is initially described 
by one of the static solutions in (14) (not necessarily the ground state solution Co), and 
in which a black hole is formed by collapsing matter fields. First we consider the simple 
shock wave solution, but our results can be easily extended to general infalling matter 
configurations. The shock wave of infalling matter is described by (T+ + ) c £ = -^rf5(x + — Xq) 
and (t£_) c £ = 0.[8] Unlike in the RST model, here we have a general initial static geometry, 
and the shock wave forms a black hole only if M, the energy of the shock wave, is above a 
certain threshold energy. We assume that M is above that threshold. Integrating (T+ + ) c £ 
in (10) and using (12) and a± = 0, we find the evaporating black hole solution 

k M 

e -2^ = e -2p = _X 2 X + X - _ _l og (_A 2 x + X-) - Yj( X+ ~ X )®( X+ ~ X o) + C ' ( 15 ) 

where Q(x) is the standard step function. 

Before the shock wave, i.e., in the region x + < Xq , we have a static solution (14) which 
is not globally flat. If |-[1 — log(«/4)] + C < 0, then the scalar curvature diverges on a 
timelike curve a = cr s , for which e" 2 ^) =0. Of course this is a region of strong coupling, 
and one would expect to have higher order quantum corrections there. On the other hand, if 
^■[1 — log(«/4)] + C > 0, the scalar curvature is bounded. Then the region {x + > 0,x~ < 0} 
is geodesically incomplete and one can analytically extend it to x~ > and x + < 0. Also 
in this case there is a region of strong coupling near a = <J m in = — 2T^°S( ^ ) • I n the 
semiclassical approximation, one avoids the strong coupling region by imposing boundary 
conditions on a suitable time-like hypersurface. [24, 14,25] For the static solutions (14), 
(Tl+ia*)) and (Tl_(a ± )) are constant on any time-like hypersurface a = const. Moreover 

{Tl+i**))^ = (TLiv^h^ for any constant <r . (16) 
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This means that we can limit our model to a region in which the semiclassical approxima- 
tion is valid by imposing reflecting boundary conditions (16) on any time-like hypersurface 
a = o"o that lies outside the region of strong coupling (these boundary conditions are 
also conformal [25]). The geometry before the shock wave is therefore a static geometry, 
defined in the region a > a s in the case f [1 — log(«/4)] + C < (or a > a m i n in the 
case — log(«/4)] + C > 0), with reflecting boundary conditions on a = cr s + 8 (or on 
c = Cmm + 8) where 8 is an arbitrary small positive constant. 

The solution to the future of the shock wave (x + > x^ ) is (see (15)) 



20 _ p -2p _ _\2 X +( X - + A) - -log(-\ 2 X + X-) + + C, (17) 



M 



e ^ = e 



X 

M 

i 

singularity at e -2 ^ = 0. The black hole singularity curve is 



where A = + . This solution is asymptotically flat and describes a black hole with a 

A x 



K M 

-X 2 xt(x: + A) - -log(-A 2 x+xJ) + — + C = 0. (18) 

Initially the singularity is behind an apparent horizon <9_|_e -2 ^ = 0,[26,[10]] which is the 
curve 

-A 2 4(x-+A) = ^. (19) 

When the apparent horizon is formed, the black hole starts radiating. One can see this by 
calculating (T^ v ) at future null infinity (x + — > oo). From (17) we see that the asymptot- 
ically flat coordinates on S" 1 " are cr 1 * 1 , related to x ^ by the conformal coordinate transfor- 
mation, Xa + = log(Xx + ) and —Act - = log( — X(x~ + A)). Using (11) and (12) we get 

ei-^n.. = ^ (i - (1 + • m 

This is the "standard" Hawking radiation in 2D, where the Hawking temperature Tjj = 
^ is a constant. [8] One can further verify that when the black hole evaporates over a 
long period of time, i.e., if M >> kX } the spectrum of the Hawking radiation is indeed 
Planckian.[[2],27] 

As the black hole evaporates by emitting Hawking radiation, the apparent horizon 
shrinks and eventually meets the singularity in a finite proper time. They intersect at (see 
Fig. 1) 

xf nf = — (e ( «a _ _) anc l x ~. = -A( 1 - -e" ( +1) ) . (21 



At this point the singularity becomes naked. We show that it is possible to impose a 
boundary condition in which a weak shock wave emanates from the intersection point, 
resulting in a solution that is stable (having non-negative ADM mass), conserves energy, 
and is continuous with the metric defined to the past of the null hypersurface x~ = x~ nt . 




Fig. 1: Penrose diagram describing formation and subsequent evaporation of a black hole in our 
model. 



Before considering the solution to the future of the null hypersurface x~ = x~ nt (the 
end-state solution), we calculate the total amount E ra( i of energy radiated during the 
evaporation. Integrating (20) over S -1 " (up to x~ nt ) gives 

(T f __(a-))dd- = M + AC - — (log(«/4) - 1) - — , (22) 

where o~ nt = &~{ x ~i n t)- The result (22) is exact. The ADM mass [22] of the dynamical 
solution (15) (relative to the ground state C = Co) is Madm = M + A(C — Co). We 
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see that the black hole radiates almost all of its initial energy. The unradiated mass 8M 
remaining as x~ — > x~ nt (which is the Bondi mass) is 

8M = M ADM - E rad = ^ (log(«/4) - 1) - \C + (23) 

We now consider the solution to the future of the point of intersection (x^~ nt , x~ nt ). A 
natural candidate for such an end-state in our model is one of the static solutions (14), 
so we try to find boundary conditions such that the solution (17) is continuously matched 
to one of the static solutions (14). Remember that the asymptotically flat coordinates are 
cr 1 * 1 , so one should replace a in (14) with a = |(cr + — &~)- In the x^ 1 coordinates the 
corresponding static solution is (see (13)) 

e -2<f> = e -2 P = _\2 X +( X - + A) - ^log(-\ 2 x + (x- + A)) + C. (24) 

We would like to see if there exists a constant C = C*, such that on the null hypersurface 
x~ = x~ nt the solutions (17) and (24) can be matched continuously. This is indeed the 
case and from (21), (17) and (24) we get C* = — -f (1 — log(«/4)). The end-state solution, 
or "remnant" , is therefore 

e -2* = e -2 P = _\2 X +( X - + A) - ^log(-\ 2 x+(x- + A)) - ^(1 - log(«/4)), (25) 
where x~ > x~ nt . From the constraint equations (9) we find that 

1 N AA 
(TL(d-)U = ? I>-/*) 2 = ~ Km)- (26) 

i = l ^ X int 

This describes a shock wave originating at the intersection point and carrying a small 
amount of negative energy, n\A/(4:X~ nt ) } to null infinity. One may call it a "thunder- 
pop". [14] The solution (25) is one of the static solutions that is asymptotically flat (with 
no radiation) on S" 1 ". This means that there is no Hawking radiation after the thunderpop 
(26). 

The mass remaining after the shockwave (26) has been emitted is 8M — n\A/(4:X~ nt ). 
One readily verifies that this is equal to the mass of the "remnant" (relative to Co) 
M rem = \(C* — Co). The fact that energy is exactly conserved, including terms of order 
h, supports the self-consistency of our semi-classical theory. Notice that C* and therefore 
the "remnant" mass is independent of the mass M of the infalling matter and of the con- 
stant C describing the initial static geometry. Even more surprising is the fact that the 
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end-state solution with C = C* is the critical solution separating singular and non-singular 
static solutions described by Eq. (24). For C > C* the curvature of the solution (24) is 
bounded, while for C < C* the curvature diverges on a time-like curve a = a S} for which 
e -2<f>{(T s ) _ q j n so l u tion space, the solution (25) that has C = C* is the boundary between 
these two different classes of solutions. 

Consider the late-time space-like hypersurface E shown in Fig. 1. Its right boundary 
(a — > oo) is i° , while its left boundary is the curve a = cr cr , for which e~ 2 ^ = 0. For 
the critical solution we have d x +(e~ 2( t >( - a ' cr ^) = e -2 ^ ") = and the curve a = o cr is 
the analytical continuation of the apparent horizon to the region x~ > x~ nt . We define 
e = a — cr cr , and calculate the metric near e = 0. From (25) we get 

2 — di 2 + de 2 
~* 2A 2 e 2 +(9(e 3 )' ( } 

where t = |(cr + + &~)- An important feature of (27) is that there is no linear term (in 
e) in its denominator. The first non- vanishing term is of order e 2 , which means that the 
geometric structure near e = is that of an infinite throat. Consider for example the 
distance along t = constant curves. The distance to a = a cr diverges logarithmically, 
exactly as it does in higher-dimensional extremal black holes. The end-state space-time 
is geodesically complete. One may consider this solution as an "extremal 2D black hole" . 
On a = a cr the Ricci scalar is constant, = 4A 2 and the geometry is regular. 

The most natural choice of Co for ground-state solution is the one with Co = C*. 
This solution describes a static radiationless geometry which is regular everywhere. Any 
solution (14) with smaller ADM mass (C < C*) has a naked singularity. In the class 
of solutions with no naked singularities, C = C* is the one with lowest energy. This is 
very similar to the linear dilaton vacuum solution (LDV) in classical dilaton gravity or to 
Minkowski space in Einstein gravity. § Also if we choose Co = C*, then the mass remaining 
after the thunderpop (26) is exactly zero. Thus the end-state solution (25) is the static 
ground-state. Its geometrical structure is independent of the initial conditions and is a 
semi-infinite throat extending into the strong coupling region. 

In our 2D semiclassical model, one does not recover all the information of the initial 



state from the end-state solution. For infalling matter described by a general (T^_ + ) c i of 



S Classical solutions with ADM mass smaller than the LDV have a naked singularity, as do 
Schwarzschild solutions with mass smaller than zero. 
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compact support, the solution (10) will depend only on the first two moments of (T|_|_) c ^, 
M = Xj x + (T^_ + ) cl dx + and P_|_ = J (T^_ + ) cl dx + .[14] The end-state solution will still be 
(25), but with A = A _2 P_|_. The information encoded in this "remnant" (or more precisely, 
in its past null boundary x~ = x~ nt ) is only about P_|_ and M. Thus in our semiclassical 
model this end-state solution does not qualify as the "cornucopion" of Ref. [28]. However, 
the semi-infinite throat extends to a region of very strong coupling. There may be sufficient 
freedom in this strong coupling region to encode more information through strong quantum 
gravitational effects. 

In this work we constructed an action in 2D dilaton gravity and showed that, with 
a natural boundary condition, all evaporating black holes in our model end in a unique 
ground-state geometry having a semi-infinite throat. 
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